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ABSTRACT 

We obtained optical spectra of the counterpart of the ultraluminous X-ray source NGC 5408 X-l 
using the FORS spectrograph on the VLT. The spectra show strong high excitation emission lines, Hen 
A4686 and [NeV] A3426, indicative of X-ray photoionization. Using the measured X-ray spectrum as 
input to a photoionization model, we calculated the relation between the Hen and X-ray luminosities 
and found that the Hen flux implies a lower bound on the X-ray luminosity of 3 x 10 39 ergs _1 . 
The [Nev] flux requires a similar X-ray luminosity. After subtraction of the nebular emission, the 
continuum appears to have a power-law form with a spectral slope of — 2.0^q'2- This is similar to low- 
mass X-ray binaries where the optical spectra are dominated by reprocessing of X-rays in the outer 
accretion disk. In one observation, the continuum, Hen A4686, and [NeV] A3426 fluxes are about 
30% lower than in the other five observations. This implies that part of the line emission originates 
within one light-day of the compact object. Fitting the optical continuum emission and archival X-ray 
data to an irradiated disk model, we find that (6.5 ± 0.7) x 10 -3 of the total bolomctric luminosity 
is thermalized in the outer accretion disk. This is consistent with values found for stellar-mass X-ray 
binaries and larger than expected in models of super-Eddington accretion flows. We find no evidence 
for absorption lines that would permit measurement of the radial velocity of the companion star. 

Subject headings: black hole physics - galaxies: individual: NGC 5408 galaxies: stellar content - 
X-rays: galaxies - X-rays: black holes 



1. INTRODUCTION 

"Ultraluminous X-ray sources" (ULXs) are variable X- 
ray sources in external galaxies with luminosities (if the 
emission is isotropic) greatly exceeding those of Galac- 
tic X-ray binaries. The irregular variability, observed 
on time scales from tenths of seconds to years, suggests 
that ULXs are binary systems containing a compact ob- 
ject. The compact objects in these binaries are either 
stellar-mass black holes in an unusual state with very 
high accretion rates or intermediate-mass black holes. It 
is of significant interest to determine the nature of the 
compact objects in ULXs. 

Beyond the luminosities, there are two lines of X- 
ray evidence suggesting that ULXs harbor intermedi- 
ate mass black holes: X-ray spectral evidence for ac- 
cretion disks co oler than those in stellar-mass black hole 
X-ray binaries ([Kaaret et al.ll2003T ) and quasiperiodic os- 
cillations (QPOs) at lower frequencies than found in 
stella r-mass black hole X-ray binaries a t s imilar luminosi- 
ties dStrohmaver fe Mushotzkvi |2003f ) . iFeng fc Kaaretl 
(|2007af) recently found that the putative cool disk in 
one ULX, NGC 1313 X-2, does not follow the rela- 
tion that the luminosity is proportional to temperature 
to the fourth power, L oc T 4 , as expected for accre- 
tion disk emission, ruling out this interpretation for that 
source. This finding casts doubt on the X-ray spectral 
evidence at least for that source. In contrast, the X- 
ray timing evidence appears robust with the QPO fre- 
quencies following the relations found in X-ray binaries 
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(|Feng fc Kaaretll2007rJ K QPOs have been definitively de- 
tecte d from only two ULXs: one in M82 and NGC 5408 
X-l (IStrohmaver fc Mushotzkvi [2001 IStrohmaver et all 
120071 ). Further, NGC 5408 X-l is s urrounded by a power- 
ful radio nebula ([Lang et al.ll2007l ) which requires an ex- 
tremely energetic outflow from the compact object. NGC 
5408 X-l is one of the best intermediate mass black hole 
candidates amongst the ULXs. 

Optical spectroscopy of ULXs provides, perhaps, the 
best means to better understand the nature of the ULXs. 
Indeed, the most definitive constraint on the compact ob- 
ject mass would come from optical measurement of the 
mass function of the companion star. The ULX in M82 
may lie in a superstar cluster and, in any event, lies in 
too crowded and obscured a region to per mit extraction 
of an optical spectrum of the ULX binary ([Kaaret et al.l 
I200ll2004al ). In contrast, NGC 5408 X-l lies in a rel- 
atively uncrowded region and is only weakly obscured. 
We recentl y identified the o ptical counterpart to NGC 
5408 X-l dLang et all [2007m . Thus, NGC 5408 X-l is 
the best intermediate mass black hole candidate ULX 
for which optical spectroscopy of the binary system can 
be performed. 

Optical emission can arise from a companion star, an 
accretion disk surrounding the compact object, or a sur- 
rounding nebula. If the light arises from the companion 
star, then one would expect continuum emission and, for 
most spectral types, a series of absorption lines. Mea- 
surement of the continuum shape would help constrain 
the spectral type of the companion star. Detection of 
spectral lines from the companion star would enable us to 
determine the star's spectral type, thereby constraining 
the evolutionary history of the binary. Repeated obser- 
vations would enable measurement of the mass function 
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Fig. 1 — Image of NGC 5408 in the F606W filter (V-band) from 
the WFPC2 on HST. The position of the optical counterpart of 
NGC 5408 X-l, the position of star A, and the FORS1 slit (0.51" 
wide) are marked. The arrow is 4" long and points north. 



TABLE 1 

VLT/FORS Observations of NGC 5408 X-l 



Observation 


Start 


Duration 


Air mass 


OBI 


2008-04-08 04:59:04 


3x850 


s 


1.050 


OB2 


2008-04-08 05:46:04 


3x850 


s 


1.049 


OB3 


2008-04-08 06:45:13 


3x850 


s 


1.095 


OB4 


2008-04-09 05:46:17 


3x850 


s 


1.051 


OB5 


2008-04-10 05:08:33 


3x850 


s 


1.046 


OB6 


2008-04-10 05:54:29 


3x850 


s 


1.056 



Note. — The table lists VLT/FORS observations of NGC 5408 and 
includes the date and UT time of the start of the observation, the 
number of integrations and the duration of each integration, and the 
average airmass during each observation. 



and a lower bound on the compact object mass. 

If the optical emission arises from the accretion disk, 
then we would expect to see a blue or hot continuum and 
emission lines. High-excitation lines, particularly Hen, 
would give direct information about the nature of the ac- 
cretion flow and the incident X-ray flux. These lines can 
provide information about the nature of the binary sys- 
tem and an estimate of the true X-ray luminosity. In ad- 
dition, since the accretion disk moves with the compact 
object, the emission line radial velocity represents the 
radial velocity of the compact object and thus provides a 
means to constrain the mass of t he compact objec t . This 
technique was demonstrated by iHutchings et al.1 (|1987l ) 
for the persistent black hole candidate LMC X-l. 

If the optical light arises from a surrounding nebula 
then the nebula could be either shock-powered, in which 
case low excitation lines would be expected, or X-ray 
photoionized, in which case high excitation lines would 
be expected. Forbidden lines can produced in a photoion- 
ized nebula but not in the high densities of an accretion 
disk. Nebulae may also produce continuum emission. 

We obtained several observations of NGC 5408 X-l 
using the FORS instrument on the VLT and analyzed 
archival X-ray data from XMM-Newton. The observa- 
tions and data reduction are described in §2. The results 
are described in §3 and discussed in §4. 



2. OBSERVATIONS AND ANALYSIS 
2.1. Optical 

The Very Large Telescope (VLT) was used for the ob- 
servations because of the dimness of the optical coun- 
terpart of NGC 5408 X-l, V = 22.4, and the need to 
obtain moderately high resolution spectra. NGC 5408 
X-l transits nearly overhead at Cerro Paranal. We 
chose to use the FOca l Reducer and Spectrograph 1 
(|Appenzeller et al.lfl998l . FORS1) which has a blue sen- 
sitive CCD because we are primarily interested in lines 
that would be produced in a photoionized nebula or ac- 
cretion disk and many of the important lines for such 
systems are in the blue range. FORS1 is mounted on the 
VLT unit telescope 2 (UT2). We used the GRIS_600B 
grism which covers a wavelength range of 3350 A to 
6330 A with a dispersion of 0.75 A/pixel and a spec- 
tral resolution of A/AA = 780 at the central wavelength. 
We used a 0.51" slit and 2x2 binning of the pixels with 
high gain readout and standard resolution. 

The spectrograph slit was oriented so that the clos- 
est neighbors of the ULX counterpart are best separated 
along the slit. The slit also covered a b right star, 'star 
A', that appears in the 2MASS catalog (|Skrutskie et al.1 
12006ft at a position of R.A. = 14 h 03 m 18'?97, Deck = 
-41°22'56'/6 (J2000). Since the ULX optical counter- 
part is quite dim, star A was used for acquisition. We 
checked the slit position by correcting the astrometry 
of the FORS acquisition image using the stars in the 
2MASS catalog. Fig . [1] shows an HST image of the field 
(|Kaaret et al.ll2004bl ) with the reconstructed slit position 
superimpos ed. The opti c al cou nterpart of NGC 5408 X-l 
identified in lLang et al.l (|2007f ) appears to be well placed 
in the slit. 

Multiple observations of NGC 5408 X-l were obtained 
over three successive nights, see TableHJ The data reduc- 
tion was performed u sing the Ima ge Reduction and Anal- 
ysis Facility (IRAF) (|Todvlll993h . Starting with the raw 
data, we created bias and lamp flat-field images using ex- 
posures obtained on the 3 nights of observations. These 
corrections were applied to the spectrum images. Each 
observation block (OB) consisted of three integrations 
of 850 s each with an offset of 12 pixels along the spa- 
tial axis between successive integrations. We analyzed 
the data by combining the 3 two-dimensional spectra in 
each observation block in order to remove bad pixels and 
cosmic ray hits. We used the xregister task in IRAF to 
caculate the offsets needed to align the images and then 
the imcombine task to average the images, applying the 
cedclip method to remove cosmic rays and bad pixels. 
We also combined the entire set of 18 spectra into one 
using the same procedure. 

The dimness of the continuum emission from the opti- 
cal counterpart of the ULX made finding a trace of the 
spectrum unreliable. Instead, we first traced the spec- 
trum of star A and then used that as a reference for the 
trace of the ULX spectrum. The trace position on the 
spatial axis varied by less than 0.5 pixels along the length 
of the dispersion axis. The trace for the ULX counterpart 
was centered using the profile of the Hen A4686 emission 
line. Wavelength lamp and standard star (EG247) spec- 
tra were used to wavelength and flux calibrate the spec- 
tra. Extinction corrections were applied using the CTIO 
extinction tables available in IRAF. 
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The spectra were searched for emission and absorption 
lines using the splot task in IRAF. A number of emission 
lines, discussed below were found. No stellar absorption 
lines were detected. There are some weak absorption fea- 
tures that are due to subtraction of night sky background 
lines at 5580 A and just above the He I line at 5885 A. 
The latter night sky line makes measurement of the He 
I A5876 flux unreliable. After correction for the redshift 
of NGC 5408, z = 0.00168785, the lines are all within 
1 A of the expected values except for [Nev] A3426 which 
differs by 3 A. However, [Nev] is at a shorter wavelength 
than the shortest available in the wavelength calibration 
lamp spectra, 3650 A, and we believe that the offset is 
due to extrapolation of the calibration. 

2.2. X-ray 

To determine the spectrum of X-rays from the ULX, 
we extracted the longest publicly available observation of 
NGC 5408 from the XMM-Newton archive. The observa- 
tion was p erformed on Jan 13, 200 6 and had a duration 
of 133 ks (|Strohmaver et all I2007T ) . We performed the 
standard XMM data filtering, including removal of data 
from times when background flares were present using 
SAS version 8.0.0. We selected only single pixel events 
in order to obtain the best energy resolution, since ad- 
equate photon statistics were available, and binned the 
X-ray spectrum in the 0.3-9.55 keV band to have at least 
25 counts per bin. We fitted the spectrum using XSPEC 
version 12.5.CK 

Following iKaaret et alj (|2004bl), we chose to use a 
Comp tonization model, compps (|Poutanen fe Svenssonl 
Il996f ). including multicolor blackbody emission from an 
accretion disk because we must extrapolate the X-ray 
spectrum to energies below the minimum observed X- 
ray energy in order to provide a useful input to the 
photoionization modeling described below. We used the 
tbabs model to account for interstellar absorption along 
the line of sight. We found evidence for an emission 
line or the combination of a unresolved group of emis- 
sion lines near 0.9 keV and included a Gaussian to fit 
that component. Our best fitted model gave a \ 2 °f 
243.2 for 219 degrees of freedom with a column density 
N H = (7.07 ± 0.03) x 10 20 cm~ 2 , an inner disk temper- 
ature kT = 172 ± 8eV, and a Comptonization electron 
temperature of kT = 74 ± 8 keV and optical depth of 
0.43 ± 0.07. We note that the electron temperature and 
optical depth depend on the specific geometry and form 
of the electron distribution used. The Gaussian emission 
had a centroid of 0.91 ± 0.02 keV, a width of 67 ± 20 eV, 
and an equivalent width of 31 cV. This emission com- 
ponent may represent a t hermal plasma, as previ ously 
suggested for this source (IStrohmaver et all I2007T) an d 
also seen from some other ULXs (|Feng fc Kaaretl 120051 ). 
However, using the apec thermal plasma model instead 
of the single Gaussian produces a worse fit, \ 2 of 263.3 for 
219 degrees of freedom. Thus, interpretation of the line 
is not clear and higher resolution spectra will likely be re- 
quired to understand the origin of the emission. Correct- 
ing for absorption and integrating the model from 1 eV 
to 100 keV, the total luminosity is 1.0 x 10 40 ergs _1 . 

We note that our X-ray and optical observations are 
not simultaneous; they are separated by two years. How- 
ever, NGC 5408 X-l does not appear to be highly vari- 
able in X-rays. The X-ray light curve constructed by 
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Fig. 2. — Spatial profiles of emission from the ULX. The solid 
curve is continuum emission in the band from 4370 A to 4680 A, the 
peak near pixel 92 is due to the optical counterpart of the ULX 
and the peak near pixel 99 is due to another star. The dashed 
curve is the Hen emission line. The dotted curve is the continuum 
emission from star A, shifted to align with the ULX profile and 
shown for comparison. The continuum emission from the ULX 
appears point like, while the Hen is slightly extended. The dot- 
dashed curve is H/3. The dot-dot-dot-dashed curve is [Oin]. Both 
of these emission lines are broader than and displaced from the 
ULX continuum. Near the bottom of the plot, the upper solid 
horizontal line shows the spatial region used for the narrow trace 
and the lower horizontal line shows the wide trace. These are the 
spatial regions used for extraction of spectra as discussed in the 
text. Each pixel is 0.25". 



TABLE 2 

Optical emission lines from NGC 5408 X-l 
Identification Wavelength (A) Flux FWHM (A) 



[Ne V] A3426 


3428.9 


315 ± 21 


5.2 


± 


0.2 


[O II] A3727 


3734.4 


1103 ± 40 


6.2 


± 


0.1 


Hz, A3771 


3777.0 


46 ± 11 


4.5 


± 


1.1 


H6» A3798 


3804.8 


56 ± 12 


5.4 


± 


1.1 


H?j A3835 


3842.4 


71 ± 11 


5.1 


± 


0.8 


[Ne III] A3869 


3875.9 


285 ± 15 


4.9 


± 


0.2 


HC, He I A3889 


3896.1 


188 ± 14 


5.5 


± 


0.3 


He, [Ne III] A3970 


3976.6 


260 ± 16 


5.8 


± 


0.3 


H<5 A4102 


4109.5 


268 ± 14 


5.4 


± 


0.2 


H7 A4340 


4348.6 


471 ± 15 


5.1 


± 


0.1 


[O III] A4363 


4371.2 


81 ± 10 


5.3 


± 


0.7 


He II A4686 


4693.7 


382 ± 12 


6.3 


± 


0.2 


H/3 A4861 


4869.9 


1000 ± 8 


5.2 


± 


0.1 


[O III] A4959 


4967.8 


1083 ± 5 


5.1 


± 


0.1 


[O III] A5007 


5015.7 


3226 ± 4 


5.1 


± 


0.1 


He I A5876 


5885.3 


57 ± 6 


4.9 


± 


0.8 


[O I] A6300 


6310.5 


101 ± 5 


5.2 


± 


0.5 



Note. — The table gives the dereddened fluxes of emission lines 
relative to the flux of H/3 = 1000. The dereddened H/3 flux is (8.6 ± 
0.9) X 10~ 16 erg cm -2 s -1 . 



IKaaret et all ([20031 ) spanning 22 years with multiple ob- 
servatories, showed a ratio between the maximum and 
minimum observed fluxes of 1.4. Those fluxes ranged 
from 7% below to 30% above our measured flux. Five 
Chandra observations carried out over 3 years showed a 
variability of 12% and no significant changes in spectral 
shape. 
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Fig. 3. — Spectrum of the optical counterpart of NGC 5408 X-l. The dereddened flux is plotted versus redshift-corrected wavelength. 



3. RESULTS 

Figure [2] shows profiles along the spatial direction of 
the H/3 emission line, the Hell emission line, and inte- 
grated flux over the 4370-4680 A continuum band. The 
profiles were extracted from the combination of the entire 
set of 18 spectra. The spatial profile of the continuum 
flux from Star A in the same band, shifted to match the 
profile of the ULX, is shown for comparison. The con- 
tinuum emission from the ULX appears unresolved. The 
Hell emission is well centered on the continuum emis- 
sion from the ULX, but is somewhat broader. We es- 
timate a spatial extent for the Hem region of roughly 
1.3" (full width at half maximum) corresponding to a di- 
ameter of 30 pc for a distance of 4.8 Mpc to NGC 5408 
(jKarachentsev et alj|2003) . The H/3 emission is offset 
from the continuum and much broader. We chose to ex- 
tract spectra from a trace that is 1.75" wide covering 
the continuum emission and the central part of the Hen 
emission. The trace is 7 pixels wide as marked by the 
horizontal line on Fig. [2l For comparison, we also exam- 
ined the spectrum from a 4.5" wide trace. The H/3 fluxes 
increases by a factor of 2.0 and the Hen A4686 increases 
by a factor of 1.3 for the wide trace. Since the nebula 
appears to have a larger angular size than our slit size 
(0.51") and the background region selected (due to the 
presence of star A) may contain some nebular emission, 
the fluxes quoted here may underestimate the true fluxes. 

Using the Balmer decrements of H7/H/3 and H<5/H/3 
to estimate the reddening, we find E(B-V) = 0.08±0.03. 
The reddening estimated from dust within the Milky 
Way is E(B - V) = 0.068 (jSchlegel. Finkbeiner. fc Davisl 
(19981) and we take this as a lower limit on the redden- 
ing. The column density measured via the X-ray spec- 
trum is equivalent to an opti cal reddening E(B-V) = 
0.13 (jPredehl fc Schmittl Il995l) . If the optical extinc- 



TABLE 3 
Emission line and continuum fluxes 



Line 


OBI 


OB2 


OB3 


OB4 


OB5 


OB6 


[O III] A5007 


223 


213 


206 


219 


228 


216 


H/3 A4861 


76 


70 


70 


73 


75 


72 


He II A4686 


25 


25 


21 


17 


26 


25 


[Ne V] A3426 


18 


19 


18 


11 


18 


20 


Continuum 


1.9 


2.0 


1.8 


1.4 


2.0 


2.1 



Note. — The tabic lists the observed fluxes of four selected emission 
lines in units of 10 -17 ergem -2 s~ 1 and the continuum flux integrated 
over the band 4000-6000 A in units of 10 — 14 erg cm -2 s~ 1 for each 
observation block (OB). 

tion measures to the edge of a nebula surrounding the 
ULX and the additional extinction seen in the X-rays is 
due to absorption in the nebula, then the column density 
through the nebula would be Ng = 2.4 x 10 20 cm~ 2 . We 
corrected for reddening using t he ex tinction curve from 
ICardelli. Clavton. and Mathisl (|l989f l with R v = 3.1. 

The combined spectrum from all of the observations 
is shown in Fig. [3] The spectrum is dominated by 
emission lines, but continuum emission is also detected. 
Table [2] gives the dereddened line fluxes relative to 
H/3 and their widths. The dereddened H/3 flux was 
(8.6 ±0.9) x 10- 16 ergcm- 2 s _1 . The widths of the lines 
are all consistent with the instrumental resolution except 
for [O II] AA3726, 3729 which is a blend, He II A4686, and 
possibly [Ne V] A3426. 

We analyzed each observation block individually. Most 
of the line fluxes appear constant with variations on the 
order of 5% that may be due to differences in the place- 
ment of the slit. However, the fluxes for He II A4686 and 
[Ne V] A3426 in OB 4 are only 70% of the average in the 
other observations, see Table [3] The average width of 
the [O III] A5007 line is 5.lA. The average width of the 
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Fig. 4. — Continuum spectra of the optical counterpart of NGC 
5408 X-l. The dereddened flux is plotted versus rcdshift-corrected 
wavelength. The crosses show the average spectrum, the diamonds 
are for OB6 that had the highest flux, and the squares are for OB4 
that had the lowest flux. The solid line is a power-law fitted to the 
average spectrum in 4000-6000A range. The data lie above the fit 
at wavelengths shortward of the Balmer edge, indicating a Balmer 
excess that may be due to nebular emission. 

He II A4686 line, excluding OB4, is somewhat broader at 
6.2A. This suggests an intrinsic width of the He II line 
of 3.5A, equivalent to 230 km/s. The width of the He II 
A4686 line in OB4 is 5.6A, which is marginally narrower 
than in the in other observations. 

We measured the continuum emission by removing all 
the detected lines and then averaging the remaining flux 
m bins of 50A. We fit a power-law to the dereddened 
and redshift-corrected continuum spectrum in the 4000- 
6000A band. The total flux in this band for each OB 
is given in Table [3J The level of continuum emission 
is roughly constant, within about 5%, except for OB4 
which has about 70% of the average flux of the other ob- 
servations. The best fitted spectral index for the average 
of all observations is — 2.01q'3,, where the uncertainty al- 
lows for the uncertainty in the reddening. The spectral 
indexes for the individual OBs are consistent with this 
range. Fig. [5] shows the continuum emission for the aver- 
age spectrum and the OBs with the highest (OB6) and 
lowest (OB4) continuum flux. 

4. DISCUSSION 
4.1. Emission lines 

The fact that the line emission appears spatially ex- 
tended and the narrow widths of the lines suggest that 
they are produced in a nebula. However, the variability 
of the He II A4686 and [Nc V] A3426 lines suggests that 
part of the emission in these two lines is produced close 
to the compact object, within one light day. 

The spectral resolution is not adequate to resolve the 
[O II] A3729/A3726 lines and the coverage does not ex- 
tend to the [S 11] A6716/A6731 lines, so these standard 
density indicators cannot be used to estimate the den- 
sity. From the flux and spatial extent of the H/3 line, a 
rough estimate of the average density is rin ~ 7cm~ 3 . 
The column density from the ULX to us through (half 



of) the nebula is then in reasonable agreement with the 
excess absorption in the X-ray spectrum. 

The weakness of low-ionization lines, particularly [O 
1] A6300 suggests that shocks are not the dominant 
source of ionization. Using the ratio of [O ill] lines, 
(/a4959 + /a5007)//a4363 = 53±7, we estimate the nebular 
temperature to be near 17000 K. This is consistent with 
that expected for photoioni zation and lower than that ex - 
pected for shock ionization (jOsterbrock fc Ferlandl2006h . 
The metallicity indicator ([O n](A3727 + [O m](A4959 + 
A5007))/H/3 = 5.4±0.2, suggesting a metallicity some- 
what lower than solar (jKewlev fc Dopital 120021 ) . Us- 
ing the spectrum from the 4.5" trace, the same ratio 
is 5.7 ±0.1. 

The strength of the Hen A4686 and [Nev] A3426 emis- 
sion is strong evidence that the line emission is pow- 
ered by photoionization from a source of high exci- 
tation. The ratio Hcii/H/3 equals 0.38 for the aver- 
age of all observations and 0.27 for OB4. Both val- 
ues are higher than found in H11 regions in nearby 
galaxies (|Garnett et al.l |1991[ ) and higher than in the 
X- ray photoionized nebul a surrounding Holmberg II X- 

I (jPakull fc Mirionil[2003l : iKaaret et al.ll2004rJ ). Within 
the 1.75" trace in the average spectrum, the luminosity 
in Hen A4686 is 9 x 10 35 ergs" 1 for a dis tance of 4.8 Mpc 
to NGC 5408 (jKarachentsev et al.ll2002] ). Taking the full 
extent of the nebula (within the slit) increases this to 
1.2 x 10 36 ergs _1 . The [Nev] A3426 luminosity integrated 
over the full extent of the nebula is 1.0 x 10 36 crgs~ 1 . 
These values may underestimate the true luminosities, 
since the nebula appears to be more spatially extended 
than our 0.51" slit. The variability of the Hen line in- 
dicates that part of the emission is due to the compact 
object instead of the extended nebula. The Hen lumi- 
nosity for the wide trace in OB4 is 1.0 x 10 36 ergs -1 . 

The Hen emission from NGC 5408 X-l is a factor of 10 
more luminous than that from the p ersistent stellar-mass 
black hole X-ray binary LMC X-l (jPakull fc AngebaultJ 
Il986t ) but similar in luminosity and size to the Hem re- 
gion ar ound the ultraluminq us X-ray source Holmberg 

II X-l (jKaaret et al.l I2004H ). There are powerful ra- 
dio nebulae associated with both ULXs (jMiller et al.l 
120051: lLang et al.l I2007T). while LMC X- l is radio quiet 
(|Fender. Southwell, and Tzioumisl Il998f ). The higher 
Hell luminosities of the ULXs suggest that the compact 
objects produce a greater isotropic luminosity than does 
the stellar-mass black hole X-ray binary LMC X-l. 

We used the photo ionization code Cloudy version 07.02 
iFerland et al.l (|1998l ) to model the nebula with the X-ray 
spectrum described above as input in order to determine 
the relation between the Hen luminosity and the ionizing 
X-ray luminosity. We placed the X-ray source at the 
center of a spherical nebula with a filling factor of 1 and 
a radius of 20 pc. The radius is slightly larger than the 
observed size of the Hem region since the H/3 emission 
extends to lower excitation and thus larger radii. 

The relation between Hen luminosity and the total lu- 
minosity of the ionizing X-ray/UV radiation is not sensi- 
tive to metallicity or density as long as the density is high 
enough that the region of fully ionized He is contained 
within the simulated nebula and the density is not so high 
that there is strong absorption of Hen A4686 within the 
nebula. From the measured Hen luminosity for OB4, the 
minimum value we observe, we estimated that the true 



6 



Kaaret & Corbel 



X-ray luminosity must be at least 2.5 x 10 39 erg s _1 . This 
value is a lower bound since an increased X-ray luminos- 
ity would be required if nebula surrounding the ULX 
covers only part of the emitted X-rays or if the radia- 
tion field at the outer edge of the nebula is still strong 
enough to ionize He. Given the X-ray flux and spectrum 
measured with XMM-Ncwton, the total intrinsic X-ray 
luminosity is 1.0 x 10 40 ergs _1 and the X-ray luminosity 
in the observed 0.3-10 keV band is 6.6 x 10 39 ergs~ 1 af- 
ter correction for absorption. Thus, the X-ray emission 
from the ULX is not strongly beamed. Use of the [Nev] 
luminosity produces similar conclusions, but they are a 
bit less certain due to the influence of metallicity. 

A Cloudy simulation with a density of njj = 6cm~ 3 , 
an input X-ray luminosity of 2.5 x 10 ergs , and a 
metallicity of 0.5Z Q produced results that matched the 
data for OB4 reasonably well. In particular, the [O in] 
line ratio, (/\4959 + /a5007)//m363 = 54. The diameter of 
the Hem region was close to 30 pc. The Hen luminosity 
was 1.0 x 10 36 ergs _1 and the [Nev] luminosity was 8 x 
10 35 ergs" 1 . All of these are close to the observed values 
for OB4. The ratio Heii/H/3 = 0.52 is somewhat higher 
than observed because our simulated nebula docs not 
cover the full Hn region. 

We note that the observed H/3 and [O in] emission is 
offset from the continuum emission, while the Hen emis- 
sion is reasonably well centered on the continuum. Thus, 
while our Cloudy simulation is likely a reasonable approx- 
imation to the geometry of the high excitation regions 
of the nebula, it does not adequately describe the more 
extended, low excitation parts of the nebula. Therefore, 
while the simulation results regarding the Hen and [Nev] 
fluxes arc likely robust, the results regarding the low ex- 
citation lines may be dependent on the geometry of the 
nebula. High resolution imaging of the nebula in several 
emission lines will be necessary to understand its struc- 
ture in detail. 

The nebula is similar in size to the radio nebula near 
NGC 5408 X-l (jLang et al.ll200l . The ULX Holmberg 
II X-l also has both Hell and ra dio nebulae and t heir 
morphology appears to be linked (|Miller et alj 120051 ). It 
would be of interest to image both the optical and radio 
in more detail and compare their morphologies. 

4.2. Continuum emission 

Compared with a power-law fitted to the data in the 
4000-6000 A band, the continuum spectrum, see Fig. [4] 
shows a marked increase at short wavelengths with higher 
flux at wavelengths shorter than the Balmer edge. Such 
a feature is expected from a photoionized nebula. The 
solid line in Fig. [5] shows the predicted continuum emis- 
sion from the Cloudy simulation described below and 
the triangles show the data after the nebular continuum 
is subtracted. The nebular continuum matches the en- 
hancement across the Balmer edge reasonably well. The 
excess in the subtracted flux just longward of the Balmer 
edge is likely due to high order lines in the Balmer series 
that the line removal procedure applied to the optical 
data does not successfully identify and remove. 

We fitted the continuum after nebular subtraction with 
a power-law in the 4000-6000 A band. The best fitted 
spectral index is — 2.0lig is an d the dominant errors are 
the nebular subtraction and the reddening correction. 
The fitted flux is 8.1 x 10~ 18 crgc m - 2 s" 1 A" 1 at 4300 A 
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Fig. 5. — Continuum spectrum of the optical counterpart of 
NGC 5408 X-l. The dereddened flux is plotted versus redshift- 
corrected wavelength. The crosses show the average continuum 
spectrum. The dashed line shows the continuum emission from the 
photoionized nebula as calculated via Cloudy. The triangles show 
the average continuum spectrum after subtraction of the nebular 
component. The solid line is a power-law fitted to the nebula- 
subtracted continuum spectrum in 4000-6000A range. The dotted 
line shows the disk irradiation model fitted to the X-ray and optical 
data. 

and 4.9 x 10~ 18 photons cm" 2 s" 1 A" 1 at 5500 A. These 
translate to Johnson-Cousins magnitudes of Bq = 22.3 
and Vo = 22.2. The latter value is in reasonable agree- 
ment with that found from HST images of the counter- 
part, Vo = 22.2 (jLang et alJl2007t ). 

The X-ray to optical flux ratio, £ = B Q + 2.5 log -Fx 
where Fx is the X-ray flux density in /iJy after correc- 
tion for absorption, is £ = 19.6 at 2 keV and £ = 22.5 
at 0.5 keV. As noted above, our X-ray and optical ob- 
servations are separated by 2 years. However, since the 
observed X-ray variability of NGC 5408 X-l is relatively 
modest, sec fl2.2| the induced uncertainty in £ is only 
0.3. These values of £ are similar to those found for 
those low-mass X-ray binaries where the optical spec- 
tra are dominat ed by reprocessing of X-rays in th e outer 
accretion disk (Ivan Paradiis fc Mc C lint o eld fl995l) . The 
spectral shape, a power-law with a spectral index is close 
to —2, is also similar to X-ray binaries where the optical 
light is primarily due to reprocessed X-rays. Thus, a sig- 
nificant fraction of the light from the optical companion 
may arise from reprocessing in the disk. 

We fitted the X-ray and optical data simul- 
taneously with the ir radia ted disk model of 
iGierliriski. Done, fc Pagel (|2008D . diskir. We first 
used the Cloudy simulation based on the X-ray spec- 
trum modeled with compps to estimate the nebular 
optical continuum. We then fitted the X-ray and optical 
nebula-subtracted continuum emission data with the 
diskir model. We then ran Cloudy again, using this 
initial set of diskir model parameters. We used this 
second Cloudy run to, again, estimate the nebular 
optical continuum; this is the dashed line in Fig. [5l 
We then performed the nebula subtraction and fitted 
the X-ray and optical data with the diskir model a 
second time. The parameters from the second iteration 



Optical light from NGC 5408 X-l 



7 



of Cloudy are very similar to those of the first iteration. 
The ratio of Hen luminosity to bolometric luminosity of 
the compact object agrees within 5% and the other line 
fluxes agree within 20%. 

Our best fitted irradiated disk model gave a \ 2 of 278.7 
for 258 degrees of freedom with a column density Nh = 
(6.35± 0.07) x 10 20 cm~ 2 , an inner disk temperature kT = 
158±8eV, an asymptotic photon index of 2.36±0.13, an 
electron temperature of kT = 1.67^5' 3± keV, and a ratio 
of Compton to disk luminosity of 0.66 ± 0.15. The frac- 
tion of the Compton tail thermalized in the disk was fixed 
at 0.1. The irradiated portion of the disk had an outer 
radius of 8600 ± 1700, in units of the inner disk radius, 
r m = 6GM/c 2 , where M is the compact object mass. 
The irradiated disk thermalized a fraction of the bolo- 
metric luminosity / ou t = (6-5 ± 0.7) x 10 -3 . Correcting 
for absorption and integrating the model from 0.1 eV to 
100 keV, the total luminosity is 9.4 x 10 39 crgs _1 . As for 
the Comptonized disk model, the irradiated disk model 
required addition of a Gaussian near 0.92 keV to produce 
an adequate fit. We note that the X-ray and optical ob- 
servations were not simultaneous. We redid the fits with 
the optical fluxes reduced by 30% and the only signifi- 
cant change in the parameters was a reduction of / out to 
4.9 x 10~ 3 . 

The fact that the optical continuum is well described 
by a model of the irradiated outer parts of a standard 
geometrically thin, optically thick accretion disk raises 
confidence in the application of standard accretion disk 
models to the spectra of NGC 5408 X-l. The frac- 
tion of the total luminosity thermalized in the outer 
disk is higher than seen from stellar-mass X-ray bina- 
ries in the soft or 'thermal dominant' spectral state, but 
is within in the range seen in states where a sizable frac- 
tion of t he emission arises from the Co mptonized com- 
ponent (|Gierlihski. Done, fc Pagal2008l ). The spectral 
state of NGC 5408 X-l would be classified as the 'steep 
power-law state' on the basis of the energy spectrum 
a nd the presence of QPOs in the classification scheme 
of iRemillard fc McClintockl (|2006h . It would be of inter- 
est to study reprocessing in stellar-mass X-ray binaries in 
similar states, in which a moderately soft power-law com- 
ponent supplies significant flux, to allow a direct compar- 
ison. From the best fit value found for the outer radius 
of the irradiated disk, we calculate an orbital velocity at 
the outer edge of the disk of 1200 km/s. This is larger 
than the measured width of the Hen emission line, which 
must arise at least in part from reprocessing since it is 
variable. These results may suggest that we view the 
disk close to face-on. 

As noted in the introduction, the cool temperatures 
derived from standard disk models applied to ULX spec- 
tra have been int erpreted as evidenc e for intermediate 
mass black holes ([Kaaret et al.1 12003T ) . but the validity 
of such models have b een ruled out for some sources 
iFeng fc Kaaretl (|2007al ). The alternative is that the 



sources are stellar mass black holes with supcr-Eddington 
accretion rates. In the so-called 'slim disk' models of 
super-Eddington accretion flows, the disk surface is con- 
vex causing the outer disk to b e shielded from the lum i- 
nous portions of the inner disk (|Abramowicz et aTlll988h . 
Thus, the fraction of the apparent luminosity that ir- 
radiates the outer disk is much lower than for a stan- 
dard thin disk. Super-Eddington accretion may, instead, 
cause the ejection of mass leading to a wind that re- 
processes the hard spectrum of radiatio n from the inner 
part of the disk into a softer spectrum ([Poutanen et al.l 
l2007f l. The outer parts of the disk would then see this 
softer spectrum, that we measure to have a characteris- 
tic temperature near 0.16 keV. Such soft X-rays do not 
thcrmalize in the disk and, thus, do not contribution sig- 
nifican tly to the reprocessed emission i n the 4000-6000A 
band (|Gierlihski. Done, fc P age 2008). Such accretion 
flows would thcrmalize only a very small fraction of the 
apparent luminosity in the outer disk. 

The measured ratio of optical to X-ray flux for NGC 
5408 X-l is larger than would be expected in models of 
super-Eddington accretion flows, but is in the range ex- 
pected from a standard thin disk. Further observations 
of the continuum emission from NGC 5408 X-l would 
allow us to better constrain the properties of the irradi- 
ated disk. Simultaneous measurement of the optical and 
X-ray spectrum would provide a better constraint on the 
thermalized radiation fraction. Extension of the optical 
coverage into the IR would enable measurement of the 
roll over of the irradiated spectrum to a Rayleigh-Jeans 
spectrum. From our best fitted irradiated disk model, 
this should occur at roughly 1 /im. Extension of the op- 
tical coverage into the UV could enable measurement of 
the long wavelength tail of the direct disk emission. 
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